November 26, 2007 


VARIATIONS OF THE MID-IR AROMATIC FEATURES 
INSIDE AND AMONG GALAXIES 

F. Galliano 

Observational Cosmology Lab., Code 665, NASA Goddard Space Flight Center, Greenbelt 

MD 20910, USA 

Department of Astronomy, University of Maryland, College Park, MD 20742, USA 

gallianoQastro . umd . edu 
S. C. Madden 

Service d’Astrophysique, L’Orme des Merisiers, CEA/Saclay, 91191 Gif-sur-Yvette, France 

A. G. G. M. Tielens 

NASA Ames Research Center, Mail Stop 245-3, Moffett Field, CA 94035, USA 

E. Peeters 

NASA Ames Research Center, Mail Stop 245-6, Moffett Field, CA 94035, USA 
SETI Intstitute, 515 N. Whisman Rd, Mountain View, CA 94043, USA 
and Physics & Astronomy Dept., University of Western Ontario, PAB 213, London ON 

N6A 3K7, Canada 

and 

A. P. Jones 

Institut d’Astrophysique Spatiale, Universite de Paris XI, 91405 Orsay, France 

ABSTRACT 

We present the results of a systematic study of mid-IR spectra of Galactic 
regions, Magellanic Hu regions, and galaxies of various types (dwarf, spiral, 
starburst), observed by the satellites ISO and Spitzer. We study the relative 
variations of the 6.2, 7.7, 8.6 and 11.3 /xm features inside spatially resolved objects 
(such as M 82, M 51, 30 Doradus, M 17 and the Orion bar), as well as among 
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90 integrated spectra of 50 objects. Our main results are that the 6.2, 7.7 and 
8.6 pm bands are essentially tied together, while the ratios between these bands 
and the 11.3 pm band varies by one order of magnitude. This implies that the 
properties of the PAHs are remarkably universal throughout our sample, and that 
the relative variations of the band ratios are mainly controled by the fraction of 
ionized PAHs. In particular, we show that we can rule out the modification of 
the PAH size distribution as an explanation of these variations. Using a few well- 
studied Galactic regions (including the spectral image of the Orion bar), we give 
an empirical relation between the I 6 . 2 /I 11.3 ratio and the ionization/recombination 
ratio Go/n e x \/T gas . Finally, we discuss the physical interpretation of the I 6 . 2 /I 11.3 
ratio, on galactic size scales. 

Subject headings: dust - HII regions - ISM: structure - galaxies: dwarf, starburst 
- infrared: general 


1. INTRODUCTION 

The reprocessing of stellar light by dust in the infrared (IR) is widely used to probe em- 
bedded star formation. However, in the absence of other constraints, the physical properties 
which can usually be derived from an almost-featureless grain continuum emission is limited 
to global quantities, such as the dust mass and its average temperature. The ubiquity of 
numerous mid-IR aromatic features, in a wide variety of astrophysical objects and environ- 
ments, potentially provides more articulate diagnostics of the physical conditions. Indeed, 
these features dominate the mid-IR spectra of evolved stars (e.g. Blommaert et al. 2005; 
Kraemer et al. 2006), the cool ISM (e.g. Abergel et al. 2005; Flagey et al. 2006; Povich et al. 
2007), as well as whole galaxies (e.g. Verma et al. 2005; Smith et al. 2007), that have been 
extensively observed by the Infrared Space Observatory (ISO), and are currently investigated 
with a higher sensitivity by the Spitzer Space Telescope. In our Galaxy, one third of the stel- 
lar light is reprocessed by dust, while this fraction can go up to 99%, in starburst galaxies. At 
solar metallicity, roughly 15% of the cooling is radiated through the most powerful mid-IR 
bands, centered at 3.3, 6.2, 7.7, 11.3 and 12.7 pm. 

Historically, these emission features were attributed to very small grains (~ 10 A), tran- 
siently heated by single photon absorption, in order to account for the independence of the 
color temperature with the distance from the illuminating star, in several reflection nebu- 
lae (Sellgren 1984). In parallel, the central wavelengths of these bands were recognized to 
coincide with the vibrational modes of aromatic material (Duley Sz Williams 1981). These 
features are now commonly attributed to the molecular modes of Polycyclic Aromatic Hy- 



-3- 


drocarbons (hereafter PAHs; Leger & Puget 1984; Allamandola et al. 1985, 1989), which 
are planar molecules made of 10 to ~ 1000 carbon atoms, excited primarily by ultraviolet 
(UV) photons. With silicate and carbon grains, PAHs are a main component of dust models 
(Desert et al. 1990; Dwek et al. 1997; Draine & Li 2001; Zubko et al. 2004). Their absorption 
efficiency has been modeled using astrophysical observations, laboratory measurements and 
quantum theory (in particular Desert et al. 1990; Joblin et al. 1992; Verstraete et al. 2001; Li 
& Draine 2001; Mattioda et al. 2005a, b; Draine & Li 2007; Malloci et al. 2007). In addition 
to being major radiative coolants of the interstellar medium (ISM), PAHs are responsible for 
most of the photoelectric heating of the gas in photodissociation regions (hereafter PDRs) 
and the neutral interstellar medium, due to their high cumulated surface area (e.g. Bakes 
& Tielens 1994; Hollenbach & Tielens 1997). For the same reason, they probably play an 
important role in the grain surface chemistry (e.g. Tielens & Allamandola 1987). In our 
Galaxy, they contain ~ 15 — 20% of the depleted carbon (Zubko et al. 2004, with solar 
abundance constraints). As a consequence, they are believed to be part of the interstellar 
carbon condensation chain (Cherchneff et al. 2000; Dartois et al. 2005). 

From an extragalactic point of view, the luminosity of the 6.2 /im feature can be used 
as a tracer of star formation (Peeters et al. 2004b). However, this tracer is biased by global 
parameters such as the ISM metallicity. Indeed, PAHs are underabundant in low-metallicity 
galaxies (e.g. Galliano et al. 2003, 2005, 2007; Draine et al. 2007). There is a general correla- 
tion between the PAH-to-continuum intensity ratio and the ISM metallicity (Madden et al. 
2006; Wu et al. 2006; O’Halloran et al. 2006; Smith et al. 2007), and consequently between 
the IRAC 8 M m/MIPS 24 M m broadband ratio and the metallicity (Engelbracht et al. 2005). The 
origin of this trend has been attributed to radiative and mechanical destruction mechanisms 
by Madden et al. (2006) and O’Halloran et al. (2006) respectively. Conversely, from the 
detailed modeling of the spectral energy distribution (SED) of nearby galaxies, Galliano 
et al. (2007) showed that the PAH-to-gas mass ratio at different metallicities coincides with 
the relative amount of carbon dust condensed in the envelopes of low-mass stars, during the 
Asymptotic Giant Branch phase (AGB). This study suggests that PAHs are injected into the 
ISM by their progenitors, the AGB stars, several hundreds of Myr after the beginning of the 
star formation, when the gas has already been enriched by more massive stars. This delay 
corresponds to the time needed for AGB stars to evolve off the main sequence. Therefore, the 
delayed injection of AGB-condensed carbon dust into the ISM offers a natural explanation 
for the paucity of PAHs in low-metallicity environments. 

From a cosmological point of view, the large luminosity in these IR emission features 
coupled with the high sensitivity of Spitzer has allowed the detection of PAHs in distant 
luminous infrared galaxies out to redshift 2 ~ 2 (e.g. Elbaz et al. 2005; Yan et al. 2005; 
Houck et al. 2005). Hence, understanding what controls the properties of the aromatic 
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features on large scales is required, in order to properly interpret broadband surveys. 

The detailed characteristics of the mid-IR features, such as their shape, their central 
wavelength or the intensity ratio between the different bands are known to vary (see Peeters 
et al. 2004a). Such variations are essentially due to modifications of the molecular structure 
of the PAHs, in different astrophysical environments. In particular, since each feature is 
attributed to a given vibrational mode, the ratio between these features will vary with 
quantities such as the charge, the hydrogenation or the size and shape of the molecule. The 

3.3 yum PAH band arises from the radiative relaxation of CH stretching modes, while the 

11.3 /m i feature originates in the CH out-of-plane bending modes; CC stretching modes are 
responsible for the features between 6 and 9 ^m; CH in plane bending excitation produces 
part of the 8.6 /im band. Now, laboratory studies and quantum calculations shows that the 
CC modes are instrinsically weak in neutral PAHs, and become stronger when the PAHs are 
ionized (Langhoff 1996; Allamandola et al. 1999; Bauschlicher 2002; Kim & Saykally 2002). 
Therefore, the 6 to 9 yum bands will be much more intense for a PAH + than for a PAH 0 , while 
it will be the opposite for the 3.3 and 11.3 yum. Consequently, the ratios between the CC 
and the CH feature intensities depend on the charge of the PAHs, which is directly related 
to the physical conditions (e.g. intensity of the ionizing radiation field, electron density, etc.) 
in the environment where the emission is originating. 

Evidence of variations between features in different astrophysical environments have 
been reported by many authors. For example, Joblin et al. (1996) showed that the I 8 . 6 /I 11.3 
ratio decreases with increasing distance from the exciting star, in the reflection nebulae 
NGC 1333 - where R is the integrated intensity of the feature centered at A [im. Hony 
et al. (2001) found a good correlation between the 3.3 and 11.3 yum CH bands, in a sample 
of Galactic Hll regions, YSOs, and evolved stars, while they reported variations of I6.2/I11.3 
by a factor of 5. The observations of Galactic and Magellanic Hu regions, presented by 
Vermeij et al. (2002), indicate that the ratios l6.2/In.3> I7.7/I11.3 and Is.6/Iii.3 are correlated. 
Furthermore, they suggest a segregation between the values of these ratios in the Milky 
Way and those in the Magellanic Clouds. Bregman & Temi (2005) studied the variation of 
I7.7/I11.3 in three reflection nebulae. Assuming that this variation is controlled by the charge 
of the PAHs, they could relate this band ratio to the ratio Go/n e between the integrated 
intensity of the UV field, Go, and the electron density, n e . Similarly, Compiegne et al. (2007), 
studying the detailed variations of the mid-IR spectrum in the Horsehead nebula, attributed 
the high relative strenght of the In .3 feature to a high fraction of neutral PAHs, due to the 
high ambient electron density. On the contrary, Smith et al. (2007) studied the variation 
of I7.7/I11.3 coming from the nuclear regions of the SINGS legacy program galaxies. They 
find that this ratio is relatively constant among pure starbursts, but varies by a factor of 
5 among galaxies having a weak AGN. They interpret this effect as a selective destruction 
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of the smallest PAHs by the hard radiation arising from the accretion disk, ruling out the 
explanation in terms of ionization of the molecules, in these particular environments. This 
interpretation is also supported by the high 11.3 jum and the weakness of the 3.3 /im band 
in the AKARI spectrum of the giant elliptical galaxy NGC 1316 (Kaneda et al. 2007). 

The previous considerations stress the diversity of the possible interpretation of the mid- 
IR feature variations in galaxies. We need to identify the main physical processes controling 
the PAH bands, if we are to use them as diagnostic tools. This is the aim of this paper. It 
presents a quantitative analysis of mid-IR spectra (ISO and Spitzer) of Galactic regions, low- 
metallicity dwarf galaxies, quiescent spirals, starbursts and AGNs. To achieve our goal, we 
focus on identifying the main trends between the various PAH features, at different spatial 
scales and in different environments. We then link these variations to the physical conditions 
inside the studied region. Preliminary results of this study were published in Galliano (2004, 
2007). 

The paper is organized as follows. In §2, we present our sample and the data reduction. 
We discuss the spectral decomposition of the mid-IR spectrum, in §3. The results of this 
decomposition are presented in §4; we study the various trends between the band ratios 
within galaxies, and among different types of environment. Then, in §5, we provide a physical 
interpretation of these trends, when the structure of the ISM is resolved, and when it is not. 
Finally, we summarize our conclusions in §6. 


2. OBSERVATIONS AND DATA REDUCTION 
2.1. The Sample 

In order to systematically study the properties of the mid-IR aromatic features, we 
include in our sample a wide variety of Galactic regions and galaxies, covering a large range of 
metallicities and star formation activities. We merge the ISO samples of starbursts and AGNs 
presented by Laurent et al. (2000), spirals by Roussel et al. (2001), Magellanic regions and 
dwarfs by Madden et al. (2006), and complement them by low-metallicity sources observed 
with Spitzer. We add to this sample, the ISO / CAM spectra of several Galactic regions, like 
M 17 (Cesarsky et al. 1996b), NGC 2023 (Abergel et al. 2002), NGC 7027 (Persi et al. 1999) 
and the Orion bar (Cesarsky et al. 2000), and the ISO/SWS spectra (2.5 — 45 pm) of several 
compact Hll regions published by Peeters et al. (2002b). The SMC Bl#l spectrum is the 
one presented by Reach et al. (2000). 

The global properties of the selected sources are presented in Table 1. If relevant, the 
distances were homogenised to H 0 — 71 h~ 2 kms -1 Mpc -1 . 
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Table 1. General properties of the sample. 


Name 

R.A. 

(J2000) 

Dec. 

(J2000) 

Aperture 

Mid-IR ] 

spectrograph 

Distance 

(Mpc) 

12 + log(0/H) 
[ref.] 

Category 

Haro 11 

00 ,l 36 m 52?5 

-33°33'19" 

10" 

Spitzer/IRS 

92 

7.9 

[2] 

Dwarf 

SMC Bl#l 

00' l 45 m 33?0 

-73°18'46" 

84" - 36" 

750/CAM 

0.06 

8.0 

[3] 

Magellanic 

NGC 253 

00 h 47 m 32^9 

— 25°17 / 18" 

40" 

750/CAM 

3.3 

9.0 

[4] 

SB/AGN 

NGC 253 p 

00 h 47 m 32!9 

-25°17'18" 

10" 

750/CAM 

3.3 

9.0 

[4] 

SB/AGN 

NGC 253 e 

00 h 47 m 32?9 

— 25°17'18" 

40" - 10" 

750/CAM 

3.3 

9.0 

[4] 

SB/AGN 

SMC N 66 

00' l 59 m 02!0 

-72°10'36" 

120" 

750/CAM 

0.06 

8.0 

[3] 

Magellanic 

NGC 520 

0l' l 24 m 34?9 

+03°47 , 31" 

30" 

750/CAM 

27 



Dwarf 

NGC 613 

0l' l 34 m 17!5 

— 29°24'58 ,/ 

60" 

750/CAM 

19 

9.2 

[5] 

Spiral 

NGC 613 p 

01 ,l 34 m 17!5 

-29°24'58" 

20" 

750/CAM 

19 

9.2 

[5] 

Spiral 

NGC 891 

02 /l 22 TO 33®4 

+42°20'57" 

200" 

750/CAM 

9.6 

8.9 

[6] 

Spiral 

NGC 1068 

02' l 42 m 40?6 

-oo°oo'47" 

60" 

750/CAM 

15 

9.0 

[7] 

SB/AGN 

NGC 1068 p 

02 ,l 42 m 40?6 

-00°00'47" 

20" 

750/CAM 

15 

9.0 

[7] 

SB/AGN 

NGC 1068 e 

02 ,l 42 m 40?6 

— 00°00'47 / ' 

60" - 20" 

750/CAM 

15 

9.0 

[7] 

SB/AGN 

NGC 1097 

02 ,l 46 m 19?l 

-30 o 16'28" 

100" 

750/CAM 

12 

9.0 

[8] 

Spiral 

NGC 1097 p 

02 ,l 46 m 19n 

-30°16'28" 

40" 

750/CAM 

12 

9.0 

[8] 

Spiral 

NGC 1140 

02 h 54 m 33®5 

— 10°01'44" 

20" 

750/CAM 

25 

8.0 

[9] 

Dwarf 

NGC 1365 

03 ,l 33 m 35?6 

-36°08'23" 

100" 

750/CAM 

19 

9.1 

[10] 

Spiral 

NGC 1365 p 

03 h 33 m 35?6 

— 36°08'23" 

40" 

750/CAM 

19 

9.1 

[10] 

Spiral 

NGC 1365 e 

03' l 33 m 35*6 


100" - 40" 

750/CAM 

19 

9.1 

[10] 

Spiral 

IC 342 

03' l 46 m 49 s 7 

+68°05 , 45" 

40" 

750/CAM 

3.8 

8.9 

[11] 

SB/AGN 

IC 342 p 

03 h 46 m 49?7 

+68°05 , 45" 

12" 

750/CAM 

3.8 

8.9 

[11] 

SB/AGN 

IC 342 e 

03 h 46 m 49?7 

+68°05 , 45" 

40" - 12" 

750/CAM 

3.8 

8.9 

[11] 

SB/AGN 

IC 342 map 

03 h 46 m 49?7 

+68°05 , 45" 

3" 

750/CAM 

3.8 

8.9 

[11] 

SB/AGN 

NGC 1569 

04 ,l 30 m 49?l 

+64°50'52" 

120" 

750/CAM 

2.2 

8.2 

[12] 

Dwarf 

NGC 1569 e 

04' l 30 m 49 s l 

+64°50 , 52" 

120" - 12" 

750/CAM 

2.2 

8.2 

[12] 

Dwarf 

NGC 1808 

05'*07 m 42?3 

— 37°30 , 47" 

50" 

750/CAM 

11 

9.1 

[13] 

SB/AGN 

Orion bar D8 

05 h 35 m 18?2 

— 05°24'40" 

14" x 20" 

750/SWS 

475 pc 



Hu region 

Orion bar D5 


-05 o 25'10" 

14" x 20" 

750/SWS 

475 pc 



H II region 

Orion bar D2 

05 h 35 m 21?4 


14" x 20" 

750/SWS 

475 pc 



H II region 

Orion bar 

5 h 35 m 20?0 


6" 

750/CAM 

475 pc 



PDR 

30 Doradus 

05 ,l 38 m 34?0 


120" 

750/CAM 

0.05 

8.4 

[3] 

Magellanic 

30 Doradus p 

05' l 38 m 34!0 

-69°05'57" 

40" 

750/CAM 

0.05 

8.4 

[3] 

Magellanic 

30 Doradus e 

05 ,l 38 m 34!0 

-69°05'57" 

120" - 40" 

750/CAM 

0.05 

8.4 

[3] 

Magellanic 

30 Doradus map 05 /l 38 m 34!0 

— 69°05'57" 

6" 

750/CAM 

0.05 

8.4 

[3] 

Magellanic 

NGC 2023 

05 h 41 m 38?3 

-02°16 , 33" 

14" x 20" 

750/SWS 

475 pc 



PDR 

II Zw 40 

05 h 55 m 42?7 

+03°23'29" 

24" 

750/CAM 

10 

8.1 

[14] 

Dwarf 
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Table 1 — Continued 


Name 

R.A. 

(J2000) 

Dec. 

(J2000) 

Aperture 

Mid-IR 

spectrograph 

Distance 

(Mpc) 

12 + log(0/H) 
[ref.] 

Category 

He 2-10 

08 h 36 m i5?2 

— 26°24'34" 

10" 

Spitzer/ IRS 

8.7 

8.9 

[12] 

Dwarf 

M 82 

09 h 55 m 51!8 

+69°40'46" 

90" 

ISO/ CAM 

3.6 

9.0 

[15] 

SB/AGN 

M 82 p 

09 /l 55 m 51®8 

+69°40'46" 

24" 

ISO/ CAM 

3.6 

9.0 

[15] 

SB/AGN 

M 82 e 

09 /l 55 m 51?8 

+69°40'46" 

90" - 24" 

ISO/CAM 

3.6 

9.0 

[15] 

SB/AGN 

M 82 map 

09 fc 55 TO 51!8 

+69°40'46" 

3" 

ISO/CAM 

3.6 

9.0 

[15] 

SB/AGN 

NGC 3256 

10 /l 27 m 5Hl 

— 43°54'17" 

24" 

ISO /CAM 

37 

8.9 

[16] 

SB/AGN 

NGC 3256 p 

10 ,l 27 m 5m 

—43° 54' 17" 

10" 

ISO/CAM 

37 

8.9 

[16] 

SB/AGN 

NGC 3256 e 

10 h 27 m 5in 

— 43°54'17" 

24" - 10" 

ISO/ CAM 

37 

8.9 

[16] 

SB/AGN 

Mrk 33 

10 /l 32 m 31®9 

+54 o 24'04" 

10" 

Spitzer /IKS 

20 

8.4 

[16] 

Dwarf 

Arp 299 

ll /l 28 m 31?0 

+58°33'39" 

24" 

ISO/ CAM 

41 



SB/AGN 

UM 448 

ll h 42 m 12?4 

+00°20'03" 

10" 

Spitzer /IKS 

70 

8.0 

[17] 

Dwarf 

IR 12331 

12 h 36 TO 01!9 

— 61°51'04" 

14" x 20" 

ISO/ SWS 

4.5 kpc 



H II region 

NGC 4945 

13 /l 05 m 26®2 

— 49°28'15" 

60" 

/50/CAM 

3.9 



SB/AGN 

NGC 4945 p 

13 fc 05 m 2®2 

— 49°28'15" 

20" 

/50/CAM 

3.9 



SB/AGN 

NGC 4945 e 

13 h 05 m 26?2 

— 49°28'15" 

60" - 20" 

/50/CAM 

3.9 

... 


SB/AGN 

Centaurus A 

13 /l 25 m 28?0 

-43°01'06" 

30" 

/50/CAM 

3.8 

~9 

[18] 

SB/AGN 

Centaurus A e 

13 h 25 m 28!0 

— 43°01'06" 

10" 

/50/CAM 

3.8 

~9 

[18] 

SB/AGN 

M 51 

13 ,l 29 m 52?7 

+47°11'43" 

140" 

/50/CAM 

8.4 

8.7 

[19] 

Spiral 

M 51 p 

13 h 29 m 52!7 

+47°11'43" 

60" 

ISO/ CAM 

8.4 

8.7 

[19] 

Spiral 

M 51 e 

13 /l 29 m 52?7 

+47°11'43" 

140" - 60" 

/50/CAM 

8.4 

8.7 

[19] 

Spiral 

M 51 map 

13 h 29 m 52!7 

+47°11'43" 

6" 

ISO/ CAM 

8.4 

8.7 

[19] 

Spiral 

M 83 

13'*37 m 00?7 

— 29°51'58" 

200" 

/50/CAM 

4.5 

9.2 

[20] 

Spiral 

M 83 p 

13 /l 37 m 00!7 

— 29°51'58" 

50" 

ISO/ CAM 

4.5 

9.2 

[20] 

Spiral 

M 83 e 

13 /l 37 m 00?7 

— 29°51'58" 

200" - 50" 

ISO/ CAM 

4.5 

9.2 

[20] 

Spiral 

M 83 map 

13 h 37 m 00?7 

— 29°51'58" 

6" 

/50/CAM 

4.5 

9.2 

[20] 

Spiral 

Circinus 

14 h 13 m 09?6 

— 65°20'21" 

40" 

/50/CAM 

4.0 



SB/AGN 

Circinus e 

14 h 13 m 09!6 

— 65°20'21" 

40" - 10" 

ISO/ CAM 

4.0 



SB/AGN 

Arp 220 

15 h 34 m 57!2 

+23°30'11" 

20" 

ISO/ CAM 

73 



ULIRG 

Arp 220 p 

15 ,l 34 m 57?2 

+23°30'11" 

10" 

ISO/ CAM 

73 



ULIRG 

IR 15384 

15 h 42 m 17!l 

— 53°58'31" 

14" x 20" 

ISO/ SWS 

2.7 kpc 



H II region 

NGC 6240 

16 h 52 m 58!8 

+02°24'06" 

20" 

ISO/ CAM 

98 



SB/AGN 

NGC 6240 p 

16 /l 52 m 58®8 

+02°24'06" 

10" 

ISO/ CAM 

98 



SB/AGN 

M 17 map 

18 h 20 m 22!0 

— 16°12'40" 

6" 

/50/CAM 

1.5 kpc 



PDR 

IR 18317 

00 

CO 

to 

CO 

— 07°54'47" 

14" x 20" 

ISO/ SWS 

4.9 kpc 



Hu region 

NGC 6946 

20 /l 34 m 51?2 

+60°09'17" 

140" 

ISO/ CAM 

5.5 

9.1 

[21] 

Spiral 

NGC 6946 p 

20 /l 34 m 51?2 

+60°09'17" 

40" 

ISO/ CAM 

5.5 

9.1 

[21] 

Spiral 
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Table 1 — Continued 


Name 

R.A. 

Dec. 

Aperture Mid-IR 

Distance 12H-log(0/H) Category 


(J2000) 

(J2000) 

spectrograph 

(Mpc) [ref.] 


/50/CAM 

750/SWS 


5.5 
700 pc 


9.1 


[ 21 ] 


Spiral 

PDR 


NGC 6946 e 

20 /l 34 m 51!2 

+60°09'17" 

140" - 40' 

NGC 7027 

21 h 07 m 01?7 

+42°14'09" 

14" x 20' 

IR 22308 

22' l 32 m 45!9 

+58°28'21" 

t— 1 

X 

to 

O 

IR 23030 

23' l 05 m 10!6 

+60°14'41" 

H-* *• 

X 

to 

0 

IR 23133 

23' l 15 m 31?4 

+61°07'08" 

0 

CM 

X 

H 

IR 23128 

23'*15 m 46?0 

— 59°03'17" 

20' 

IR 23128 p 

23' l 15 m 46?0 

-59°03'17" 

10' 

NGC 7714 

23' l 36 m 14?l 

+02°09'19" 

10' 


750/SWS 

5.5 kpc 

Hu region 

750/SWS 

5.2 kpc 

H 11 region 

750/SWS 

5.5 kpc 

H 11 region 

750/CAM 

180 

SB/AGN 

750/CAM 

180 ••• 

SB/AGN 

Spitzer / IRS 

37 8.5 

[22] Dwarf 


References. — [2] Bergvall et al. (2000); [3] Dufour et al. (1982); [4] Zaritsky et al. (1994); [5] Alloin et al. (1979); 
[6] Otte et al. (2001); [7] Dutil & Roy (1999); [8] Storchi-Bergmann et al. (1995); [9] Heckman et al. (1998); [10] Roy 
& Walsh (1997); [11] Pilyugin et al. (2004); [12] Kobulnicky & Skillman (1997); [13] Ravindranath & Prabhu (2001); 
[14] Perez-Montero & Diaz (2003); [15] Boselli et al. (2002); [16] Mas-Hesse & Kunth (1999); [17] Izotov & Thuan 
(1998); [18] Schaerer et al. (2000); [19] Bresolin et al. (2004); [20] Webster & Smith (1983); [21] Kobulnicky & 
Johnson (1999); [22] Gonzalez-Delgado et al. (1995). 


Note. — The sources are ordered according to their right ascension. For ISO / CAM sources, the apertures are 
circular, centered on the coordinates listed in the second and third columns. For ISO / SWS sources, the beam is 
14" x 20" (bands 1 and 2). The diameter of the aperture is given in the fourth column. When the diameter has 
the form 8 i" — Q^', it means that the we subtracted to the flux in the circular aperture of diameter #i", the flux in 
a concentric circular aperture of diameter 62 " < #i". The sources which are followed by “map” are those for which 
the signal-to-noise of the spectral-image is good; in this case, the aperture refers to the pixel field of view. The 
letter after the name of the source designates the type of aperture: when there is no letter, most of the emission of 
the object is encompassed in the aperture; “p” is for peak emission; “e” for extended emission. 
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2.2. ISO/ CAM Data Reduction 

Most of the sources in Table 1 were observed with ISO / CAM (Cesarsky et al. 1996a) on 
board the ISO satellite (Kessler et al. 1996). We refer to the work of Madden et al. (2006) 
for a detailed description of the data reduction applied to all of the ISO / CAM data used 
in this study. The CVF performed spectral imaging using a 32 x 32 detector array, with a 
sampling of 3"pixel -1 or 6 /; pixel -1 in our cases, from A = 5 p m to 16.5 fxm with one pointing 
of two CVFs, from A = 5 to 9.5 pm and from A = 9.0 to 16.5 pm. The spectral resolution 
increases from A/AA = 35 to 51 across the full spectra. 

For the data treatment, we used the CAM Interactive Reduction (CIR, version AUG01; 
Chanial 2003). The subtraction of the dark currents was performed using the Biviano et al. 
(1998) model which predicts the time evolution for each row of the detector, taking into 
account drifts along each orbit and each revolution. We masked the glitches using multi- 
resolution median filtering (Starck et al. 1999) on each block of data after slicing the cube. 
Additional deglitching was performed manually, examining the temporal cut for each pixel. 
We corrected the systematic memory effects using the Fouks-Schubert method (Coulais & 
Abergel 2000). We computed a hybrid flat-field image placing a mask on the source and 
computing a flat field outside this mask from the median of the temporal cut for each pixel. 
For the pixels which were on-source, the flat-field response was set to the corresponding cal- 
ibration flat- field. The conversion from Analog to Digital Units to mJy/pixel was performed 
using the standard in-flight calibration data base. To remove the sky contribution, sources 
smaller than the array were masked and, for a given wavelength, the median of the pixels 
which are off-source were subtracted from each pixel. For the more extended sources, we 
subtracted an independently observed zodiacal spectrum. The contribution of this spectrum 
was a free parameter varied in order to match the properly sky subtracted fluxes in the LW2 
(6.7 /im) and LW3 (14.3 pi n) broadbands. The final product is a 3D spectral-image. We 
integrated the spectrum into an aperture encompassing the entire galaxy, and obtained a 
ID spectrum. When the angular size of the source was larger than the array, we scaled the 
spectrum to match the IRAS^m broadband flux. 


2.3. Spitzer/ IRS Spectrum Extraction 

Several of the low-metallicity sources in Table 1 were not observed by ISO/ CAM, thus we 
complemented our database with mid-IR spectra from the Spitzer/IRS spectrometer (Houck 
et al. 2004b) on board the Spitzer Space Telescope (Werner et al. 2004), when these data 
were released. Among these galaxies, the spectrum of SBS 0335-052 has been published 
by Houck et al. (2004a), NGC 7714 by Brandi et al. (2004), NGC 5253 by Beirao et al. 
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(2006), Haro 11 and I Zw 18 by Wu et al. (2006). We refer to these studies for a detailed 
presentation of their mid-IR properties. We considered only low-resolution data, taken with 
the SL (Short-Low) module, from A = 5.2 pm to 14.5 pm, and the LL (Long-Low) module, 
from A = 14.0 pm to 38.0 pm, both with a spectral resolution of A/AA ~ 64 — 128. 

We retrieved the Basic Calibrated Data (BCD) that have been preprocessed by the 
Spitzer Science Center (SSC) data reduction pipeline. They have been converted to flux 
density, corrected for stray light and flatfielded. The extraction of the spectra from the 2D 
space/wavelength images was performed with the Spectral Modeling, Analysis and Reduction 
Tool (SMART, version 5.5.6; Higdon et al. 2004). We first inspected the BCD images and 
identified the hot pixels which had not been masked out by the SSC. We replaced them 
by the median of their neighbors. For each module, an off-source position is observed. We 
subtracted this spectrum from that of the source, in order to remove the sky emission. 
The extraction of the ID spectrum was performed inside a column whose width vary with 
wavelength. We have excluded the bonus order. Then, the various frames, for each nod 
position, were coadded. Since the long wavelength end of the SL module and the short 
wavelength end of the LL module were not systematically overlapping, we finally scaled the 
SL module, in order to obtain a continuous spectrum. This scaling factor can be as large 
as 50%. At the time of this publication, the Spitzer / IRS data handbook recommends that 
the signal to noise ratio not be derived from the uncertainties generated by the pipeline. 
Instead, we adopted the recommended systematic error of 20%. To take into account these 
statistical variations, we smoothed the spectra into a AA ~ 0.2 /im window (4 points), and 
took the standard deviation inside this window as the error. As we did with the ISO / CAM 
spectra (§2.2), in order to compensate the fact that we may be overlooking some extended 
emission, we scaled the spectrum to match the IRAS^m and IRAS^bum- In the case of 
II Zw 40, the Spitzer / IRS spectrum does not exhibit the PAH features that Madden et al. 
(2006) detected in the extended emission, since the Spitzer/IRS observations were conducted 
in staring mode. Indeed, this Spitzer/IRS observation does not encompass the region where 
the tail was detected. Thus we prefer to use the ISO/ CAM data for this galaxy. 


3. THE SPECTRAL DECOMPOSITION METHODS 

3.1. Inventory of the Physical Components 

We restrain our study to the spectral range 5 — 16 pm, since it is the overlap between 
the various data sets we have compiled. The main physical components contributing at these 
wavelengths are the followings. 
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• The ionic gas emits strong fine-structure ionic lines. The [Ne III] i5.56/xm> [Ne II] 12 . 81 ^ > 
[S iv] io. 5 i/im) [Ar in] 8 . 99 M m and [Ar ll] 6 . 98/im are the brightest lines. 

• Several ubiquitous broad features, originating from the molecular modes of stochasti- 
cally heated PAH molecules. The most prominent of them are centered at 6.2, 7.7, 8.6, 
11.3 and 12.7 fim. The [Neii]i 2 . 8 i/im fine is blended with the 12.7 yum feature. 

• Very small grains (VSG), fluctuating in temperature around a few hundred degrees 
Kelvin, produce a continuum emission. In addition, hot grains in Hll regions, at 
thermal equilibrium with the radiation field, may contribute to this continuum. 

• The extinction feature at 9.7 /zm, attributed to amorphous silicates, can be seen toward 
the most embedded sources. 

• Molecular lines are also present, especially the H 2 0-0 S(3)9.7 Mm line. However, they are 
much weaker than the other components. 

Fig. 1 shows the variations of the aromatic feature spectrum among galaxies (top panel), 
inside one galaxy (middle panel) and from a theoretical point of view (bottom panel). This 
figure shows that the 6.2, 7.7 and 8.6 /zm features are qualitatively tied together, and that 
most of the variations in the mid-IR aromatic spectrum is a variation of these three features 
relative to the 11.3 /xm band. The lower panel of Fig. 1 suggests that most of these variations 
could be explained by a variation of the neutral-to-cationic-PAH ratio. In what follows, we 
quantify these variations in our sample, in order to explore the validity of this explanation. 

3.2. Measuring the Intensity of the Aromatic Features 

Measuring the intensity of the aromatic features is uncertain, due to the intrinsic width 
(AA ~ 1 /zm) and the complexity of the band profiles. Indeed, contrary to gas lines, a 
large fraction of the energy of the PAH bands is radiated in the wings. These wings extend 
far outside of the central wavelength and can be difficult to reliably disentangle from the 
underlying continuum emission. Therefore, a proper extraction of the feature requires an 
assumed band profile. However, the actual profile of each band is not known. It is asymmetric 
and varies as a function of environment (e.g. Peeters et al. 2002a; van Diedenhoven et al. 
2004). Consequently, several profiles have been used in the literature, motivated by different 
physical arguments. For example, Boulanger et al. (1998) proposed that the width of the 
band is a consequence of the continuous redistribution of intramolecular vibrational energy, 
between different excitation levels, at high temperature. They estimated the transition 
timescales to be < 10 -13 s, and showed that the PAH bands were well represented by 



- 12 


[Aril] [Arlll] [SIV] [Nell] [Neill] 



Fig. 1. — Various spectra normalized to the 11.3 /mi feature. The top panel shows the total 
integrated spectra of a low-metallicity galaxy (UM 448), a normal starburst (M 82), and a 
LIRG (NGC 4945). The middle panel shows the spectra of three different regions inside M 82. 
The bottom panel shows, for comparison, the absorption coefficient of neutral and cationic 
PAHs measured in laboratory by Allamandola et al. (1999). The solid vertical lines mark 
the wavelengths of the brightest lines, and the dashed vertical lines mark the wavelengths of 
the major bands. In this figure and in what follows, F v { A) is the monochromatic flux density 
at wavelength A fim. 
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Lorentzian profiles. Conversely, a Drude profile was used by Li & Draine (2001). Such a 
profile describes the electric permeability in a solid or a large molecule. Finally, numerous 
studies simply measure the tip of the aromatic band, and assumes that the total flux radiated 
in the band scales with this quantity. 


Spline method 

[Aril] [Arlll] [SIV] [Nell] [Neill] 



Lorentzian method 

[Aril] [Arm] [SIV] [Nell] [Neill] 



Fig. 2. — Demonstration of the two spectral decomposition methods for the total spectrum 
of M 82 and 30 Doradus. 

In order to take into account the uncertainty of the actual shape of the aromatic bands, 
we systematically analysed our spectra using two different spectral decomposition methods. 
These two methods differ in the way the continuum and the aromatic bands are fitted. 
They are demonstrated in Fig. 2, on a PAH-rich spectrum (M 82), and PAH-poor one 
(30 Doradus). 

The Spline method (hereafter identified by the letter <S), used e.g. by Hony et al. (2001) 
and Vermeij et al. (2002), takes into account only the tip of the aromatic bands when 
computing their intensities, as shown in left panel of Fig. 2. The continuum (grey 
filled curve) is fitted with a spline function, constrained at the wavelengths A = 5.04, 
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5.47, 5.84, 9.18, 10.85, 13.82, 14.60, 15.08, 16.00, 16.15 /im. A second spline function 
is fitted to the spectrum, corresponding to the previous continuum, plus the “ plateau ” 
under the 6.2, 7.7 and 8.6 /im bands, and the plateau under the 11.3 and 12.7 /im 
bands. In addition to the previous wavelengths, this spline function is also constrained 
at A = 6.61, 6.80, 7.13, 8.24, 8.80, 11.85, 12.20, 13.22 /tm. The difference between the 
total spectrum and this continuum+plateau component defines the lines and features 
(line filled areas of the left panel of Fig. 2). This method is relatively robust. However, 
it does not allow us to separate the [Ne Il]i 2 . 8 i/jm from the 12.7 /im band (most of our 
spectra have a low spectral resolution). 

The Lorentzian method (hereafter identified by the letter £), used e.g. by Boulanger et al. 
(1998), Laurent et al. (2000) and Verstraete et al. (2001), assumes that the aromatic 
bands have lorentz profiles. The wings of these lorentzians account for the underlying 
plateau. A variation of this method, implemented by Smith et al. (2007), consists of 
modeling the PAH features with Drude profiles. We proceed as following. 

1. The [Nein]i 5 . 56 M m, [Ne 1 1 ] i 2 . 8 i/jmi [S iv] io.5i^rn> [Ar 1 1 1] 8 . 99 /j.ttt. i and [Ar n]6.98^m ionic 
lines are fitted with gaussian functions (diagonal line filling in the right panel of 
Fig. 2). For each (number i) of the Ny me lines, the central frequency, z/ bne , is fixed, 
the width is fixed by the spectral resolution of the instrument, Ai/ spectro , and the 
total flux of each line, i^ hne , is the free parameter. 

2. The PAH bands are fitted with lorentzian functions (horizontal line filling in the 
right panel of Fig. 2). The considered features are centered at: A = 5.3, 5.7, 6.2, 
7.7, 8.6, 11.3, 12.0, 12.7, 13.6 and 14.3 /tm. For each (number % ) of the iVband 
band, the central frequency, zA and , the width, Ai/ band , and the total flux, F band , 
are the free parameters. We emphasize the fact that we fit only one band around 
7.7 /tm, where there are actually two bands centered at 7.6 /tm and 7.8 /im. We 
proceed this way, in order to keep the number of parameters reasonably low. We 
will use the variation of the centroid of the 7.7 /im feature, noted A 7 . 7 , as an 
indicator of the 7.6 /tm to 7.8 /tm ratio. 

3. The VSG continuum is fitted with the sum of N cont modified black-bodies having 
the absorption efficiency, Qabs(^)> of graphite in the Rayleigh approximation (grey 
filled area on the right panel of Fig. 2; Laor & Draine 1993). The temperatures, 
T) cont , and the total fluxes, F £ ont , are the free parameters. We adopt N c ont = 2. 

4. The sum of the previous components is multiplied by exp(— r), where t{v) is the 
Mathis (1990) extinction law, with the Dudley & Wynn-Williams (1997) silicate 
features. The column density is the only free parameter. 
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All these components are fitted simultaneously as part of one mathematical function: 


F„(v) = (F““( v) + F„ b " d ( u) + FT») x exP (-rW) , (1) 
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the various F v being the monochromatic flux densities. An advantage of the Lorentzian 
method over the Spline method is that it allows us to separate the [Ne Il]i2.8i/xm line and 
the 12.7 fim feature, and to study the variations of the centroids of the bands. However, 
when the PAH-to-VSG ratio is very low, as in the case of 30 Doradus (Fig. 2), the 
width of the features is more uncertain. 


In both cases, I cont is defined as the integrated intensity of the continuum between 10 and 
16 /mi, and Ipah, the sum of the intensities of all the bands between 5 and 16 /mi. We 
emphasize that our methods work automatically, without any by eye adjustment. This 
condition prevents systematic effects that could result from arbitrary choices of parameters. 
These two methods are orthogonal but give similar trends (see §4 and Smith et al. 2007). 
The Spline method is clearly continuum-biased in its assumptions whereas the Lorentzian 
method is a line-biased assumption method. 

Our two methods can be applied to an integrated spectrum, or to each pixel of a spectro- 
image, provided that the signal-to- noise ratio is sufficient. We have systematically applied the 
two methods to all the spectra presented in Table 1. In our analysis, we will systematically 
consider the results of the two methods before drawing conclusions on the PAH properties. 
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4. SYSTEMATIC ANALYSIS OF THE MAJOR BAND RATIOS 
4.1. Correlations Exhibited Among Integrated Spectra 

Figs. 3-4 show select correlations between band ratios of the integrated spectra of our 
sample (Table 1), obtained with the two methods presented in §3.2. We focus on the four 
brigthest bands at 6.2, 7.7, 8.6 and 11.3 pm, and study the correlations between the various 
ratios. The intensities of the features are reported in App. A (Table 4), and the parameters 
of the various correlations are given in Table 2. As mentioned in §3.2, the two methods are 
uncertain for very low values of the PAH-to-VSG ratio. Thus we define two subgroups of 
data. 

1. The data that we consider to be reliable are marked with black symbols in Figs. 3-4. 
We define them as the measurements which have a signal-to-noise ratio, at A = 7.7 pm, 
larger than 6, and IpAH/Icont > 0.5 with the Lorentzian method , and IpAH/Icont > 0.3 
with the Spline method. 

2. The fits that we consider to be less certain are the complementary data (grey symbols 
in Figs. 3-4), having lower signal-to-noise and PAH-to-VSG ratios. 

Fig. 3 shows that there is an excellent linear correlation between the ratios l 6 . 2 /In. 3 > 
I 7 . 7 /I 11.3 and I 8 . 6 /Iii. 3 ) using either of the two methods. The variations of these ratios spread 
roughly over one order of magnitude, while the ratios I 6 . 2 /I 7 . 7 ? and I 7 . 7 /I 8.6 are roughly 
constant within the error bars (Fig. 4). In general, the measure of the intensity of the 

8.6 pm band is less accurate than for the other main bands. Indeed, it is less intense than 
the 6.2 and 7.7 pm features. In addition, it is merged with the long wavelength wing of the 

7.7 pm band. Finally, this particular feature is significantly affected by the silicate extinction 
feature around 9.7 pm, when the source is deeply embedded. This extinction is not corrected 
in the case of the Spline method, and is corrected very simply in the case of the Lorentzian 
method. That is the reason why correlations involving the 8.6 pm feature are always more 
dispersed than the others. Quantitatively, the correlation coefficients are around 0.8 for the 
four panels of Fig. 3, except for that involving the Is.6 feature, with the Spline method. Most 
of the outsiders are fits that we consider to be uncertain. 

Fig. 5 shows the consistency between the two methods. It demonstrates that, whatever 
method we use to measure the band ratios, the order of the various measures is conserved; 
i.e. a spectrum A having a lower I6.2/I11.3 ratio than a spectrum B, using the Spline method, 
will also have a lower ratio than B using the Lorentzian method. 
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These first relations, established on integrated spectra, indicate that the properties 
of the PAHs throughtout different types of galaxies and Galactic regions are remarkably 
homogeneous. They are consistent with a significant variation of the 6.2, 7.7 and 8.6 /xm 
features relative to the 11.3 /xm band, coupled with an absence of significant variations 
among the 6.2, 7.7 and 8.6 /xm bands. 


4.2. Spatial Variations of the Band Properties 

We now study the spatial variations of the properties of the mid-IR features, by ana- 
lyzing the spectral maps which have the highest signal-to- noise ratio. Similar to §4.1, we 
systematically compare the results of the two methods, but instead of applying them to 
integrated spectra, we fit the spectrum of each pixel of the spectral maps. We then degrade 
the images of each spectral component to the spatial resolution of the longest wavelength of 
the spectrum (FWHM ~ 9"). We achieve this by convolving the image of a given component 
(e.g. I 6 . 2 ), with the PSF at the longest wavelength (i.e. A = 16 /xm) deconvolved beforehand 
by the PSF of the component (i.e. A = 6.2 /xm). Figs. 6 to 11 (as well as Figs. 20 to 34 in 
App. A) show the same correlations as in Figs. 3-5, obtained inside resolved sources, instead 
of integrated spectra. Each data point represents the fit of the spectrum of one pixel. We 
selected the points, according to the criterion defined in §4.1: the signal-to- noise ratio at 
A = 7.7 /xm must be larger than 6, and IpAH/Icont > 0.5 with the Lorentzian method , and 
IpAH/Icont > 0.3 with the Spline method. The pixels which do not fall into this category 
have not been considered here. The parameters and statistical quantities relative to these 
correlations are given in Table 2. Figs. 12 to 13 (as well as Figs. 35 to 39 in App. A) show 
the spatial distributions of the components. 

In general, the correlations presented in Figs. 3-5 still hold inside individual objects 
as well. Therefore, the variations of the PAH band ratios are independent of the spatial 
resolution. For example, the pixel size is ~ 0.1 pc in M 17 and the Orion bar, while it is 
~ 1 kpc in external galaxies. From the images, we notice that the ratio I6.2/I11.3 (as well as 
the ratios I7.7/I11.3 and I8.6/I11.3) is roughly correlated with the intensity of the PAH emission, 
the highest ratios being found in the brightest regions (spiral arms, starburst region, etc.), 
while the lowest ratios are generally found in the extended, low-luminosity regions. We do 
not detect any significant variations of the centroid of the 7.7 /xm feature inside each source. 
However, the average centroid varies from one source to another. 

App. A comments on each individual source. Here we will illustrate the systematic 
variations within the individual sources by focusing on two extragalactic sources: the edge- 
on irregular starburst galaxy M 82 and the face-on spiral galaxy M 51. For M 82 (Fig. 12), 
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Lorentzian Method 



Fig. 3. — PAH properties for integrated spectra. The black error bars are for the mea- 
surements that we consider to be reliable, while the grey are more uncertain. In each 
panel representing a relation Y versus X , the grey filled area is the domain between 
Y = ay/x x X - 1 a Y /x, and Y = a Y / x x X + lcr Y /x, Y — a Y / x x X being the linear 
correlation of the data, and a Y /x, the dispersion of the data around this correlation. 
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Lorentzian method Spline method 



Fig. 4. — PAH properties for integrated spectra (continued). The same symbol conventions 
are adopted as in Fig. 3. 



l 

I 6 g/I U3 (Lorentzian) 


Fig. 5. — Comparison between the two methods. The same symbol conventions are adopted 
as in Fig. 3. 
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the maximum values of the I7.7/I11.3 ratio are found along the disc, in the star forming region, 
where the infrared emission is the highest. The value of the band ratio drops by a factor of 
~ 3 at ~ 200 pc above and below the disc, in the halo. For M 51 (Fig. 13), the I7.7 /I11.3 ratio 
is maximum along the circumnuclear ring. It is somewhat lower in the 10" nuclear region. 
The band ratio is also high along the spiral arms, coinciding with the infrared bright blobs, 
but can exhibit variations by a factor of ~ 2. Finally, the ratio drops in the interarm region. 

The highest values of the I7.7 /I11.3 — 7 ratio ( Lorentzian method 1), found in the nuclear 
star forming regions of M 82 and M 51, are similar to the values found in the compact 
Galactic Hu regions IR 15384 and IR 18317 (Table 4). Intermediate values of I7.7/I11.3 — 3 
(Lorentzian method), similar to the Galactic reflection nebula NGC 2023 are found in the 
halo of M 82 and along spiral arms of M 51. The lowest values of I7.7/I11.3 < 1 ( Lorentzian 
method), seen in the planetary nebula NGC 7027, correspond to the outermost regions of 
the two galaxies. 



M 82 (Lorentzian method) 


M 82 (Spline method) 


Fig. 6. — PAH band ratios within M 82. The grey filled areas are the correlations obtained 
for the integrated spectra (Fig. 3). The correlations obtained inside M 82 are represented 
by the hatched region. The white symbol is the value of the global measurement over the 
entire galaxy. 
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M 82 (Lorentzian method) M 82 (Spline method) 



Fig. 7. — PAH band ratios within M 82 (continued). 


M 82 



0.1 1 

Ig.g/Iii.3 (Lorentzian) 


Fig. 8. — Comparison between the two methods in M 82. The same symbol conventions are 
adopted as in Fig. 6. 
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M 51 (Lorentzian method) M 51 (Spline method) 



Fig. 10. — PAH band ratios within M 51 (continued). 


M 51 



0.1 l 

1 6.2/1 1 13 (Lorentzian) 


Fig. 11.— Comparison between the two methods in M 51. The same symbol conventions 
are adopted as in Fig. 6. 
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Fig. 12 . — Spatial distribution of the PAHs in M 82 . For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 

M 51: I, AHtot & l PAH (6.2)/l PtH (11.3) 


Lorentzian 


‘J2000 


J2000 


Fig. 13 . — Spatial distribution of the PAHs in M 51 . For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 




Table 2. Parameters of the various correlations. 
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Table 2 — Continued 


Source 

X 

Y 

Method 

ay/x 

by/x 

&Y/X 

Py/x 


I7.7/Ill.3 

A7. 7 \^ m ] 

c 

0 

7.71 

0.016 


30 Doradus 

l6.2/Ill.3 

l7.7/Ill.3 

s 

1.43 

0 

0.79 

0.78 

(Fig. 26) 



c 

2.31 

0 

0.67 

0.68 


l8.6/Ill.3 

l7.7/Ill.3 

s 

11.48 

0 

1.61 

0.35 




c 

2.10 

0 

0.62 

0.49 


I6.2/Ill.3 

l6.2/Ill.3 

S/C 

2.01 

0 

0.86 

0.44 


l7.7/Ill.3 

A7.7 

c 

0 

7.66 

0.016 


M 17 

I6.2/Ill.3 

l7.7/Ill.3 

s 

1.96 

0 

0.33 

0.81 

(Fig. 32) 



c 

2.41 

0 

0.45 

0.93 


l8.6/Ill.3 

l7.7/Ill.3 

s 

11.91 

0 

0.81 

0.33 




c 

4.04 

0 

0.84 

0.80 


I6.2/Ill.3 

I6.2/Ill.3 

S/C 

1.01 

0 

0.25 

0.61 


l7.7/Ill.3 

A7.7 [pm] 

c 

0 

7.69 

0.012 


Orion bar 

l6.2/Ill.3 

I 7 . 7 /I 11.3 

s 

1.77 

0 

0.32 

0.83 

(Fig. 29) 



c 

2.32 

0 

0.21 

0.95 


l8.6/Ill.3 

I 7 . 7 /I 11.3 

s 

8.38 

0 

0.38 

0.91 




c 

3.72 

0 

0.52 

0.86 


l6.2/Ill.3 

I 6 . 2 /I 11.3 

S/C 

1.01 

0 

0.29 

0.40 


I7.7/Ill.3 

A7.7 [pm] 

c 

0 

7.68 

0.011 



Note. — This table gives, for two given measured quantities X and Y, the 
parameters corresponding to the fit of Y by ay/x x X + by/x • The dispersion of 
the measures around this correlation is cry/x = yf {{X ~ a Y/x x X + by/x) 2 )- 
In the case where by/x = 0, py/x is the linear correlation coefficient. 
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5. INTERPRETATION OF THE AROMATIC FEATURE VARIATIONS 
5.1. Origin of the Band Ratio Variations 

The main conclusions of §4 are that (i) the ratios 16.2/111.3) I7.7/I11.3 and I8.6/I11.3 are 
correlated and span one order of magnitude throughout our sample, and (ii) the ratios 
I7.7/I6.2 and I8.6/I6.2 are not significantly anticorrelated with the I7.7/I11.3 ratio, and do not 
show significant variations within the sources which have the best signal-to-noise ratios 
(especially M 82). As explained in §1 and demonstrated in Fig. 1, the intensity of the 
features between 6 and 9 yum relative to the 11.3 /im band are one order of magnitude higher 
for PAH + than for PAH 0 . Therefore, the universal correlations displayed from Fig. 3 to 
Fig. 31 can be attributed to variations of the charge of the carriers of the aromatic features. 
However, before concluding, we first need to explore the ability of other physical processes 
to reproduce these trends. 

First, dehydrogenation of the PAHs have an effect similar to ionization, on the mid-IR 
spectrum, as shown by several laboratory studies (Allamandola et al. 1989; Allain et al. 1996; 
Jochims et al. 1999). However, as reported by Allain et al. (1996), only PAHs containing less 
than Nq — 50 carbon atoms can be considerably dehydrogenated. At the same time, this 
threshold corresponds to the minimum size of PAHs that can survive in most PDRs (Allain 
et al. 1996). For comparison, the PAH size distribution derived by Zubko et al. (2004), for 
the diffuse Galactic ISM, has a lower cut-off of N™ in = 20. Their emission is dominated by 
PAHs of Nq — 50 — 500, depending on their excitation rate. This conclusion is supported 
by the observed linear correlation of the 3.3 /xm band with the 11.3 /xm band in Galactic 
sources (Hony et al. 2001). If dehydrogenation were important then a non-linear relation 
would be expected as duo and trio groups (with bands longwards of 11.3 /xm) were converted 
into solo groups (with a 11.3 /xm band; Schutte et al. 1993; Hony et al. 2001). Likewise, 
if dehydrogenation were important, then the decrease in the I 11 . 2 /I 6.2 ratio (decreasing H 
coverage) would be expected to be accompanied by a decreasing I 12 . 7 /I 11.2 ratio (conversion 
of duo/trio’s into solo’s) in contrast to the observations (Hony et al. 2001). Therefore, we 
can neglect the effect of dehydrogenation on our mid-IR spectra. 

Second, deep extinction by the silicate feature at 9.7 yum can cause variations of the 
band ratios, as proposed by Peeters et al. (2002b), Spoon et al. (2002) or Brandi et al. 
(2006). Indeed, the 11.3 /xm feature is more absorbed than the 7.7 ytm band. This effect 
would also cause the 8.6 /urn feature to suffer from the same amount of extinction as the 
11.3 /xm band, and these two bands should be correlated. However, this is not what Figs. 3 
to Fig. 33 show. Moreover, the Lorentzian method corrects the bands for extinction (§3). 
The good agreement between the Lorentzian method and the Spline method tells us that 
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extinction can not explain the majority of the observed band ratio variations. Finally, there 
is no evidence for silicate absorption in the spectra of the Orion bax. Indeed, the total gas 
and dust column along the line of sight for this source is insufficient to produce noticeable 
10 [x m absorption. 



6 8 10 12 14 6 8 10 12 14 

Wavelength [jam] Wavelength [/um] 

Fig. 14. — Effects of the temperature distribution on the PAH spectra. The various curves 
are the theoretical emission spectra, varying the shape (Galactic or young starburst) and 
the intensity of the ISRF (xisrf), as well as the lower cut-off of the size distribution 
The size distribution is the BARE-GR-S of Zubko et al. (2004). The absorption efficiencies 
are taken from Draine & Li (2007). 

Third, a modification of the temperature distribution of the PAHs can affect the relative 
weight of the aromatic features. For example, truncating the highest temperatures would 
have the effect of decreasing the intensity of the short wavelength bands relative to the long 
wavelength ones. There are essentially two ways of modifying the temperature distribution 
of the molecules: (i) by varying the interstellar radiation field (ISRF) intensity or shape, 
and (ii) by changing their size distribution. The latter effect could happen in the vicinity 
of a very hard radiation source, like an AGN, where the smallest PAHs could undergo 
photosublimation, as proposed by Smith et al. (2007). 

We performed theoretical modeling of the PAH emission spectra, in order to study the 
latter effects. These spectra, which take into account the stochastic heating of the particles, 
are shown in Fig. 14. We adopted the PAH size distribution by Zubko et al. (2004, bare 
grains, solar abundance constraints; BARE-GR-S), and the absorption efficiencies by Draine 
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& Li (2007), with an ionization fraction of 50%. We consider successively two classes of 
ISRF, in order to explore the effect of the hardness of the radiation on the PAH emission: 
(i) the Galactic ISRF of Mathis et al. (1983) and (ii) a very hard ISRF, corresponding to 
an instantaneous burst of star formation with a Salpeter IMF, synthesized with the stellar 
evolutionary model PfiGASE (Fioc & Rocca-Volmerange 1997). The latter is normalised to 
the intensity of the Mathis et al. (1983) ISRF. In addition, we also vary the intensity of 
each ISRF, by multiplying them by a factor xisrf = 1 — 10 5 . Finally, we explore the effect 
of a possible small PAH segregation, by varying the value of the lower cut-off of the PAH 
size distribution. This lower cut-off is N ™ 111 = 20 carbon atoms, for the Zubko et al. (2004) 
model. We vary it up to N™ m = 10 3 . 

10 


n 

I s - 

I s -' 


l 


Fig. 15. — Theoretical modeling of the effect of the PAH temperature distribution on their 
band ratios, with a fixed ionization fraction (50%). The two polygons, on the right panel, 
are the variations of the band ratios with the lower cut-off of the PAH size distribution and 
with the intensity of the ISRF. The black polygon corresponds to the Galactic ISRF, and 
the grey one to the young star burst. The various arrows show the sense of variation of the 
band ratios with each parameter. 

Fig. 15 synthesizes the results of this modeling. The left panel of Fig. 15 demonstrates 
that effect of the radiation field intensity on the I 7 . 7 /I 11.3 band ratio, for the two types of 
radiation fields and the two extreme PAH size cut-offs. This figure shows that, for a given 
size distribution and ISRF type, the I 7 . 7 /I 11.3 ratio is independent of the intensity of the 
ISRF, up to Xisrf — 10 4 . Indeed, the PAHs are stochastically heated below Xisrf — 10 4 , 
therefore the shape of their emission spectrum is independent of the ISRF. However, when 
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Xisrf ^ 10 4 , the PAHs reach thermal equilibrium, and the ratio between the bands is 
modified. The right panel of Fig. 15 shows how the three effects considered here affect the 
relation between I7.7/I11.3 and 17 . 7 / 16 . 2 - The arrows indicate the direction toward which these 
parameters affect the relation. This figure shows that the combination of all these effects can 
not explain a variation of the I7.7/I11.3 ratio by more than a factor ~ 3. Moreover, it shows 
that if these effects were responsible for the variation of the I7.7/I11.3 ratio, then the I7.7/I6.2 
ratio would also be anticorrelated with it. The ratio I7.7/I6.2 is particularly very sensitive to 
Nq' u . According to Fig. 15, an increase of the I7.7/I11.3 ratio by a factor ~ 3 would imply a 
decrease of the I7.7/I6.2 ratio by a factor of ~ 2. 

One of the strong results of §4.1 and §4.2 is that the ratios I 7 . 7 /I 6.2 do not show any 
trend with I 7 . 7 /I 11 . 3 . Therefore, the correlations discussed in §4 can not be attributed to a 
modification of the temperature distribution. In particular, we can rule out an interpretation 
of the variation of the I 7 . 7 /I 11.3 in terms of destruction of the smallest PAHs. These results 
are in good agreement with those by Hony et al. (2001), who studied a wide sample of 
Galactic sources and showed that the I 3 . 3 /I 11.3 ratio was not varying significantly, while the 
I 6 . 2 /I 11.3 ratio was varying by a factor of 5. Consequently, our work shows that the variation 
of the PAH band ratios, throughout different environments, and at different spatial scales, is 
primarily controled by the ionization fraction of the molecules. In addition, it appears that 
the mixture of PAH molecules, in all these environments, is remarkably universal. 

We emphasize that our sample contains only a few AGNs. Therefore, our results reflect 
primarily the properties of star forming regions. Observations of AGN environments by 
Smith et al. (2007) and Kaneda et al. (2007) show PAH properties that differ from our 
trends, suggesting that the AGN could alter the PAH composition. 


5.2. Relating the Band Ratios to the Physical Conditions 


In §5.1, we showed that the band ratio variations were essentially due to a variation of the 
fraction of ionized PAHs in the beam. The mid-IR feature ratios are therefore related to the 
physical conditions where the emission is originating. Indeed, the fraction of ionized PAHs, 
in a given region, depends on the quantity Go/n e x (Tielens 2005, for a review), where 

Go, is the UV radiation field density, n e , the electron density, and ^ga sj the gas temperature. 
The UV field density is the integration between A = 0.09 /im and A = 0.2 /j,m of the 
monochromatic mean intensity, J\, and is normalised to the solar neighborhood value: 


r0.2 fim 
J 0.09 


Go 


47r J\ dA 


1.6 x 10- 6 Wm- 2 ‘ 


(3) 
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In this section, we use the most well-studied Galactic regions of our sample, in order to derive 
an empirical relation between the I 6 . 2 /I 11.3 band ratio and the quantity Go/n e x ^/T gas . 


Region 

Go T gas [K] 

n h [cm 1 * 3 

Reference 

NGC 2023 

1.5 x 10 4 

750 

10 5 * 

Steiman-Cameron et al. (1997) 

NGC 7027 

6 x 10 5 

2000 

10 7 

Justtanont et al. (1997) 

Orion bar (position 4) 

4 x 10 4 

500 

5 x 10 4 

Tauber et al. (1994) 


Table 3: Physical conditions in select regions. The last columns refers to the study where 
these quantities have been derived. 


Bregman & Temi (2005) gave such a calibration, from the spectro-images of three re- 
flection nebulae. In order to derive the physical conditions in their nebulae, they made 
several hypotheses. They (i) assumed a constant electronic density throughout their re- 
gions, (ii) they neglected dust extinction, and (in) they derived the relative geometry of the 
star and nebula from the profile of the scattered light. In our analysis, we have adopted a 
slightly different strategy. We have made the conscious decision to include only the sources 
for which the physical conditions have been reliably determined, using PDR models, con- 
strained by several gas lines. The sources of our sample which have been modeled in detail 
are NGC 2023, NGC 7027 and position 4 of the Orion bar (Roche et al. 1989). Their 
properties are summarized in Table 3. 


We adopt the values given in Table 3 and assume that the electrons result from photo- 
ionization of carbon and that all the gas phase carbon is ionized in PDRs, so that: 


7l e 





x n h — 1.6 x 10 4 Tin, 


(4) 


where we have adopted the carbon abundance by Sofia et al. (2004). For NGC 2023 and 
NGC 7027, we do not have reliable spatial information, hence we will only use the global 
values of the ratios. In the Orion bar, we study the spatial variations by proceeding as 
follows. 


1. We exclude the region north of the bar (see Fig. 38), which is the location of the Hu 

region cavity. Indeed, the PAH-to-continuum ratio is very low there, thus the fits are 

uncertain. Moreover, the spherical geometry of the region causes confusion along the 

lines of sight. To the contrary, the region south of the bar, is an edge-on PDR, where 

the conditions are believed to be roughly homogeneous along each line of sight. Indeed, 

it is observed to be well stratified in the main PDR tracers (CO data as well as H 2 , 
Cl, Cn and PAH data; Tielens et al. 1993; Tauber et al. 1994; Wyrowski et al. 1997; 

and HCN data by Fuente et al. 1996). 
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Orion Bar (Spline method) 



Orion Bar (Lorentzian method) 



Fig. 16. — Top panels: variation of the color of the VSG continuum (ratio of the flux at 8 
and 14 /im) as a function of the intensity of the VSG continuum, inside the Orion bar. These 
figures show that the color of the VSG continuum is virtually independent of its intensity, 
since these grains are stochastically heated. Bottom panels: variation of the PAH-to-VSG 
ratio as a function of the intensity. These figures show that the relative strength of the PAHs 
in the PDR decreases slightly toward the Hll region (when I cont increases), as an effect of 
their photo-depletion. 
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Fig. 17. — Theoretical modeling of the continuum intensity as a function of the UV radiation 
field, Go, in black. L cont / A^dust is the luminosity of the graphite component of Zubko et al. 
(2004), integrated between 10 and 16 //m, and normalized to the dust mass. The dotted 
line varies as Go- The dashed-dotted line is a power-law fit to the black line, in the range 

10 2 < G 0 < 10 4 . 
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2. We assume that the column density is homogeneous in the edge-on PDR, so that 
we can use the intensity of the VSG continuum, in order to scale the intensity of 
the UV radiation field, Gq. Indeed, the top panels of Fig. 16 show that the VSGs 
are stochastically heated. Therefore, their intensity scales with the radiation field 
density. Moreover, the intensity of the VSG continuum is probably a better tracer of 
the radiation density than the PAHs, since the latter are gradually depleted in direction 
of the Hu region (bottom panels of Fig. 16). Fig. 17 illustrates this fact by showing 
the theoretical variation of the integrated intensity of the VSG continuum between 10 
and 16 fj,m, as a function of Go- To perform this simulation, we used the radiation field 
defined in Eq. (3), and the Galactic graphite size distribution of Zubko et al. (2004, 
with solar abundance constraints). For Go < 10 3 , the mid-IR continuum intensity 
scales perfectly with Go, since the grains are out of thermal equilibrium. At higher Go, 
the relation between Go and 7 con t becomes non-linear due to temperature effects. In 
the range 10 2 < Go < 10 4 , which are the conditions in the Orion bar, the UV radiation 
field is very well approximated by Go oc I coif ■ This method has the advantage of 
taking into account the actual variations of Go, with both the distance from the star 
cluster, and the dust extinction. 

3. In order to improve the signal-to-noise ratio of our spectra, we average the pixels into 
30 bins of energy (/ con t)- The two top panels of Fig. 18 show the variation of the 
I 6 . 2 /I 11.3 ratio, as a function of the continuum intensity, with the two methods. 

4. Finally, we assume a constant electron density and temperature throughout the PDR. 


The two bottom panels of Fig. 18 shows the final calibration of the band ratio, with the two 
methods. As found by Bregman &; Temi (2005), the ratio I6.2/I11.3 (or I7.7/I11.3) increases in 
the range 10 2 < Go/n e < 10 3 . The grey stripes are the following linear fits: 


~ 1.09 x log ^ 


I11.3 

16.2 

I11.3 


L gas 


~ 1.85 x log 


n e [cm -3 ] V 10 3 K 
Go 


n e [cm -3 ] V 10 3 K 


L gas 


— 2.26 ± 0.10 ( Spline method), 


4.71 ± 0.24 ( Lorentzian method), 


( 5 ) 


valid in the range 400 



< 4000. 
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Orion Bar (Spline method) Orion Bar (Lorentzian method) 






Fig. 18. — Empirical calibration of the band ratio. The two upper panels show the variation 
of the band ratio inside the Orion bar, as a function of the VSG intensity. The black error 
bars correspond to the edge-on PDR, and the grey error bars, to the H II region cavity. The 
two lower panels show the final calibration, including NGC 2023 and NGC 7027. We kept 
only the edge-on PDR measurements, in the Orion bar. 
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5.3. Diagnosing Galaxies Using PAH Band Ratios 

In principle, Eq. (5) can be used to derive the ratio Go/n e x \/T gas from any mid-IR 
spectrum. However, if several regions with different physical conditions are encompassed 
within the spectrograph’s beam, the situation is more complex. The band ratio is then 
averaged over these various environments and its global value does not reflect a single Go/n e x 



Fig. 19.— Determination of Go/n e x ^/T gas for the sources in our sample. These figures are 
similar to Fig. 3, where the I 6 . 2 /I 11.3 has been converted into Go/n e x \jT gas using Eq. (5). 
The same symbol conventions are adopted as in Fig. 3. The small additional symbols are 
the spatial variations within M 82. 

In the case of the global spectrum of a star forming region, we can expect the I7.7/I11.3 
ratio to depend on the age and initial mass function of the stellar cluster, that will determine 
the Go at the edge of the PDR, as well as on the geometry of the ISM (dumpiness, densities), 
that will determine the transfer of the UV light into the PDR, and the electron density. If this 
hypothesis is correct, the value of I7.7/I11.3 can potentially be used in combination with other 
PDR tracers, like the [Cii]i 58mto , [Oi] 6 3 M m> [Si n] 34 . 82 ^m and H 2 mid-IR rotational lines, to 
determine the geometry and physical conditions of an unresolved star forming region. Fig. 19 
demonstrates the relation between the observed band ratios and the averaged Go/n e x \/T gas 
quantity. 

It has long been surmised that the IR emission features provide a clear mid-IR signature 
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of the interaction of FUV photons with cloud surfaces and hence a probe of the importance 
of (massive) star formation in a region (Genzel et al. 1998; Peeters et al. 2004b). The 
present study has extended this by developing the observed 6.2 to 11.3 pm band ratio as 
a quantitative tool to probe the physical conditions (e.g. Go/n e x in the emitting 

regions. We expect that this study will be of fundamental value for the interpretation of 
Spitzer data as well as future Sofia and JWST observations of galaxies in the nearby and 
early universe. 


6. SUMMARY AND CONCLUSION 

We presented the results of a systematic study aimed at understanding the main proper- 
ties of the mid-IR features in different astrophysical environments. It is based on observations 
of Galactic regions and galaxies of various types, observed by the satellites ISO and Spitzer. 
We have developed two distinct methods of spectral decomposition with different hypothe- 
ses, in order to test the robustness of our trends, and overcome eventual biases. These two 
methods have shown similar trends between the physical quantities that we have studied. 
We explored the variations of the different mid-IR features between the wavelengths 5 and 
16 pm, among integrated spectra of galaxies, and inside galaxies and Galactic regions. Our 
main results are the followings. 

1. We find that the 6.2, 7.7, 8.6 /am features are essentially tied together, while the ratio 
of these bands to the 11.3 pm feature can vary by one order of magnitude, in our 
sample. These variations are seen both inside individual sources (like M 82, M 17, 
the Orion bar, 30 Doradus, etc.), as well as among integrated spectra. In general, the 
I 6 . 2 /I 11.3 ratio is spatially correlated with the power radiated by the PAHs. It indicates 
that the ratio I6.2/I11.3 is higher in regions of intense star formation. 

2. With the help of a stochastic heating model and realistic absorption efficiencies, we 
show that the variations of the mid-IR features are essentially due to the fraction of 
ionized PAHs. We conclude that the properties of the PAHs, throughout our sample, 
are remarkably universal. In particular, we rule out the modification of the grain size 
distribution as an explanation of these variations. Indeed, we show that a modification 
of the size distribution could explain the observed variation of the I6.2/I11.3 ratio. How- 
ever, it would cause the I6.2/I7.7 and I6.2/I8.6 (and I3.3/I11.3) ratios to vary significantly 
and to be correlated with the I6.2/I11.3, contradicting our observations. 

3. The universality of the properties of the PAHs and the fact that the band ratios are 
mainly sensitive to the charge of the molecules allow us to use these band ratios as a 
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tracer of the physical conditions inside the emitting region. Using a few well-studied 
Galactic regions (including the spectral image of M 17), we give an empirical relation 
between I6.2/I11.3 and the ratio Go/n e x ^T gas . 

4. In the case where several regions with different physical conditions are integrated, the 
band ratios are dependent on the morphology of the ISM, as well as of the evolutionary 
properties of the illuminating stellar cluster. 


A. ADDITIONAL TABLE AND FIGURES 

Notes on Individual Sources 

M 82. (Figs. 6-7) Fig. 6 shows that the correlations inside M 82 are remarkably well ar- 
ticulated. The correlation coefficients (Table 2) are 0.97, for the correlations between 
I6.2/I11.3 and I7.7/I11.3. The band ratios inside M 82 follow the general trends observed 
for the integrated spectra very well. The centroid of the 7.7 /im feature is essentially 
fixed throughout M 82 (Table 2). The ratio I6.2/I11.3 follows the spatial distribution 
of the PAH emission very well, which is maximum in the central bar, where the star 
formation occurs. With good confidence (S/N > 6), we detect low ratios in the outer 
regions. 

IC 342. (Figs. 20-21) The morphology of the PAH emission (Fig. 35) follows the 12 CO(J=l— >0)2.6mm 
emission which traces out the nuclear ring and ridge associated with the stellar bar 
Sakamoto et al. (1999). The spectra of this galaxy are much noisier than those of M 82. 

The correlations are very well defined, with a linear correlation coefficient of ~ 0.8 (Ta- 
ble 2). The deviation of the relation involving the 8.6 /am feature with the Lorentzian 
method (middle right panel of Fig. 20), compared to the general correlation, is due to 
a systematic over-estimation of the attenuation, in the central region. The average 
centroid of the 7.7 /im feature is identical for IC 342 and M 82 ((A7.7) ~ 7.68 //m; 

Table 2). 

M 51. (Figs. 9-10) This galaxy is more quiescent than M 82 and IC 342. Its band ratios 
show less variation (Fig. 9) throughout the galaxy, although they are in good agreement 
with the general correlations of Fig. 3. The correlation coefficients are ^ 0.6 (Table 2) 
slightly lower than those for M 82 and IC 342, essentially because of the lower signal-to- 
noise ratio of the spectra. The spatial distribution of the I6.2/I11.3 ratio (Fig. 13) follows 
the PAH emission remarkably well enhanced along the spiral arms and the bulge. This 
band ratio reaches its maximum on the central ring and decreases somewhat in the 
center. 
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M 83. (Figs. 23-24) M 83 is a spiral galaxy which is similar in several aspects to M 51. Its 
band ratios (Fig. 23) and the spatial distribution of its PAH emission (Fig. 36) have 
the same properties qualitatively. The average centroid of the 7.7 pm is identical to 
that for M 51 and M 83 ((A7.7) ~ 7.72 p m; Table 2). 

30 Doradus. (Figs. 26-27) The band ratio correlations measured inside 30 Doradus are sig- 
nificantly dispersed (Fig. 26). The results of the 2 different feature extraction methods 
are not particularly consistent. This is certainly due to the very low PAH-to-continuum 
ratio. Indeed, Fig. 37 shows that the spatial distribution of the PAH emission, obtained 
with the Lorentzian method , is different from that obtained with the Spline method. In 
particular, the right panel of Fig. 37 shows 3 bright spots along the bar of 30 Doradus, 
that are not seen on the left panel. After verification, these are fits of such low PAH- 
to-continuum ratio pixels that the widths of the PAH bands are not constrained at all. 
The Lorentzian method determines very wide bands whose integrated intensities are sig- 
nificantly overestimated. For the same reason, the highest values of the I6.2/I11.3 ratio, 
with the Spline method (~ 2; Fig. 37), correspond to low values, with the Lorentzian 
method (< 0.3). 

M 17. (Figs. 32-33) We encounter the same trouble in the Hu region part of M 17 (upper 
left corner of the images in Fig. 39), as in 30 Doradus. For the same reasons, the 
I6.2/I11.3 ratio is maximum in this region with the Spline method (~ 2.5; Fig. 39), 
while it is minimum with the Lorentzian method (< 0.7). However, the PDR part has 
relatively high PAH-to-continuum ratios. The band ratio correlations are well defined 
inside M 17, except for 8.6 /am, with the Spline method (Fig. 32). 

The Orion bar. (Figs. 29-30) The PAH-to-continuum ratio is on average higher in this 
region than in M 17. Thus the fits are better, even along the ionization front. The 
correlation between band ratios (Fig. 29) are relatively good with linear correlation 
coefficients around 0.9 (Table 2). The plot of the consistency between the two methods 
(Fig. 31) shows a few outlying points. These fits are those of the pixels located inside 
the Hll region (upper part of the images in Fig. 38). The spatial distribution of the 
I6.2/I11.3 band ratio decreases on average from the north, where the ionizing stars are 
located, to the south with the two methods. The average centroid of the 7.7 pm feature 
is similar in M 17 and the Orion bar ((A7.7) ~ 7.68 pm; Table 2). 
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Note. — We consider that a given fit has a “high” reliability, if S/N > 6 at A = 7.7 ^m, and IpAH/^cont > 0.3 with the Spline method , and IpAH/^cont > 0.5 with the Lorentzian 
method . Otherwise, the reliability is “low”. 
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IC 342 (Lorentzian method) IC 342 (Spline method) 



Fig. 21. — PAH band ratios within IC 342 (continued). 


IC 342 



0.1 1 

Ie.a/Iu.3 (Lorentzian) 


Fig. 22. — Comparison between the two methods in IC 342. The same symbol conventions 
are adopted as in Fig. 6. 












6.2 


-48 


M 83 (Lorentzian method) M 83 (Spline method) 



I7.7/I11.3 I7.7/I11.3 

Fig. 24. — PAH band ratios within M 83 (continued). 


M 83 



0.1 1 

I e. 2/1 11. 3 (Lorentzian) 


Fig. 25. — Comparison between the two methods in M 83. The same symbol conventions 
are adopted as in Fig. 6. 
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30 Dor (Lorentzian method) 30 Dor (Spline method) 



Fig. 27. — PAH band ratios within 30 Doradus (continued). 


30 Dor 



0.1 1 

1 6 . 2 / 1 ii . 3 (Lorentzian) 


Fig. 28. — Comparison between the two methods in 30 Doradus. The same symbol conven- 
tions are adopted as in Fig. 6. 
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M 17 (Lorentzian method) M 17 (Spline method) 



Fig. 33. — PAH band ratios within M 17 (continued). 


M 17 



0.1 l 

I©. 2/1 1 1.3 (Lorentzian) 


Fig. 34. — Comparison between the two methods in M 17. The same symbol conventions 
are adopted as in Fig. 6. 
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30" 
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50 s 45 s 

“J2000 
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Fig. 35.— Spatial distribution of the PAHs in IC 342. For each method, the image is the 
total PAH intensity, and the contours are the I 6 . 2 /I 11.3 ratio. The shaded circle indicates the 
beam size. 
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Fig. 36. — Spatial distribution of the PAHs in M 83. For each method, the image is the 


total PAH intensity, and the contours are the 
beam size. 
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Fig. 37. — Spatial distribution of the PAHs in 30 Doradus. For each method, the image is 
the total PAH intensity, and the contours are the I 6 . 2 /I 11.3 ratio. The shaded circle indicates 
the beam size. 
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Fig. 38 . — Spatial distribution of the PAHs in Orion bar. For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 



Fig. 39 .— Spatial distribution of the PAHs in M 17 . For each method, the image is the 
total PAH intensity, and the contours are the I6.2/I11.3 ratio. The shaded circle indicates the 
beam size. 
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